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CONVERSION  FACTORS ,  METRIC  (SI)  TO  l’ .  S.  CUSTOMARY  AND  U.  S. 
CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


Units  of  measurement  used  in  this  report  can  he  converted  as  follows: 
_ Multiply _  By  _  _ _  To  Obtain _ 


Metric  (SI)  to  U.  S.  Customary 


nili micrometres 
micrometres 
mill  met  res 
met  res 

square  metres 

cubic  metres 

knots  (international) 

mi  1 1 irad ians 
radians  per  second 
Kelvins 


incites 

feet 

degrees  (angular) 


3.280839  x  10 
3.280839  x  10 
0.03937007 
3.280839 
10.78391 
6.102376  x  10 
1.151543 

0.05729578 

57.29578 

1.8 


-6 

-3 


feet 
f  eet 
inches 
feet 

square  feet 

cubic  inches 

miles  (U.  S.  statute) 
per  hour 

degrees 

degrees  per  second 
Fahrenheit  degrees* 


U .  S .  Customary  to  Metric  (SI) 

25.4 
0.3048 
0.01745329 


millimetres 

metres 

radians 


btuin  Fahrenheit  reading:  fro-  Kelvins,  u.- 


F  =  1 . 3 (K  -  273.15) 
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FEASIBILITY  OF  MONITORING  FLOW  PATTERNS  AND  SEDIMENT 

and  pollutant  dispersion  oi-  water  bodies" 

WITH  24-CHANNEL  SPECTRAL  DATA 


PART  I :  INTRODUCTION 


Background 


1.  Within  the  la9t  several  years,  new  iiationa^  priorities  have 
made  necessary  the  acquisition  and  analysis  of  environmental  data  on  a 
scale  not  contemplated  heretofore.  The  regions  being  considered  are  so 
large  and  the  environmental  data  requirements  are  so  universal  that 
conventional  data  acquisition  systems  are  no  longer  adequate  to  meet  the 
demands.  Even  systems  based  on  aerial  photography  and  manual  interpre¬ 
tation  are  inadequate  in  many  situations,  because  processing  is  too 
time-consuming  and  costly. 

2.  In  an  attempt  to  meet  this  challenge,  personnel  of  the  U.  S. 

Army  Engineer  Waterways  Experiment  Station  (WES) ,  sponsored  jointly  by 

the  National  Aeronautics  and  Space  Administration  (NASA)  and  the  Army 

Corps  of  Engineers,  developed  largely  automated  procedures  for  using  the 

digital  image  data  acquired  by  the  Earth  Resources  Technology  Satellite 

(ERTS-1)  to  obtain  certain  kinds  of  environmental  data  on  a  regional 

scale.  In  general,  the  procedures  based  on  ERTS-1  data  made  it  possible 

to  very  rapidly  produce  accurate  and  reliable  maps  of  open-water  surfaces 

and  of  the  concentrations  of  suspended  materials  in  the  surface  and 

1  2 

near-surface  waters.  ’  Since  the  analytical  procedure  for  mapping  the 
concentrations  of  suspended  materials  depended  upon  the  use  of  the 
spectral  data  defined  by  the  four  ERTS-1  spectral  bands,  it  was  believed 
that  similar  procedures  could  be  used  to  map  the  distributions  of  any 
other  features  of  the  landscape  that  exhibited  unique  spectral  "signatures. 
Such  features  night  be  expected  to  include  land  uses,  crop  inventories, 
a  limited  amount  of  plant  species  discrimination  (e.g.  forest  types, 
aquatic  plant  species,  etc.),  certain  kinds  of  topographic  expressions, 
and  so  on.  Thus,  a  general  analytical  procedure  for  exploiting 
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spectral  information  would  be  capable  of  supplying  a  significant  part  of 
the  environmental  data  required  for  Corps  planning  purposes. 

3.  Despite  the  fact  that  the  ERTS-1  data  proved  to  be  extraordi¬ 
narily  useful  in  many  contexts,  it  was  clear  from  the  beginning  thac 

there  were  situations  in  which  the  data  would  be  inadequate.  There  are 

3 

two  primary  reasons.  First,  the  ERTS  (and  LANBSAT*)  pixel  is  so  large 
(approximately  56.8  by  78.7  m**)  that  small  terrain  features  are 
obscured.  Since  many  items  of  interest  to  the  Corps  of  Engineers,  such 
as  small  changes  in  the  positions  of  shorelines,  are  much  smaller  than 
the  ERTS-1  pixel,  the  satellite  view  would  not  meet  all  Corps  data 
requirements . 

4.  Second,  the  ERTS-1  nultispectral  scanner  divides  the  visible 
and  near-infrared  spectrum  into  only  four  relatively  broad  bands  (i.e. 

0.5-0. 6,  0.6-0. 7,  0.7-0. 8,  and  0. 8-1.1  pm).  This  means  that  the  ERTS-1 
spectral  data  cannot  be  used  to  discriminate  between  two  features  that 
exhibit  spectral  reflectances  that  differ  only  within  one  band.  For 
example,  two  tree  species  might  exhibit  identical  spectra  as  sensed  by 
ERTS-1,  but  a  close  examination  of  each  individual  spectrum  might  reveal 
that  the  energy  from  species  A  was  concentrated  in  the  0.5-0.6-pra  range 
at  0.51  pm,  while  that  from  species  B  was  concentrated  at  0.56  pm. 

S5r.ce  this  kind  of  situation  is  common,  it  seemed  obvious  that  a  multispectral 
scanner  that  broke  the  visible  and  near-infrared  spectrum  into  narrower 
wavelength  bands  would  provide  additional  discriminatory  power., 

5.  In  view  of  the  factors  discussed  above,  the  conclusion  was 
reached  that  great  long-range  benefit  would  be  achieved  if  the  concepts 
developed  during  the  ERTS-based  research  (see  paragraph  2)  could  be 
modified  and  expanded  tc  use  digital  data  obtained  with  instruments  such 
as  the  aircraft-mounted,  NASA-operated,  24-cbanneI  Bendix  scanner 
system.  In  the  development  of  such  a  capability,  three  general  problem 
areas  were  anticipated. 


*  Formerly  ERTS. 

**  A  table  of  factors  for  converting  metric  (SI)  units  of  measurement 
to  U.  S.  customary  units  and  U.  S.  customary  units  to  metric  (SI) 
units  is  given  on  page  4. 
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a.  The  ERTS-1  (and  LANDSAT)  and  Bendix  24-channel  scanner* 

digital  computer-compatible  tapes  (CCT)  use  quite  different  data 
storage  format^,  and  thus  new  translation  algorithms  would 
be  required  in  order  to  read  the  24-channel  data, 
b^.  Because  of  differences  in  sensor  operating  altitudes  and 
mechanical  system  arrangements  and  scanning  geometry,  the 
distortions  inherent  in  the  recorded  data  are  of  different 
kinds  and  degrees.  Thus,  algorithms  developed  to  rectify 
ERTS-1  data  would  have  to  be  modified,  and  in  some  cases, 
entirely  new  algorithms  would  have  to  be  developed, 
c.  The  spectral  bands  recorded  by  the  BMSS  are  quite  different 
from  those  recorded  bv  the  EMSS,  and  the  calibration 
arrangements  by  which  sensor  output  voltages  can  be  trans¬ 
formed  into  radiance  values  are  quite  different  from  those 
in  the  EMSS;  thus,  new  analytical  algorithms  would  have  to 
be  written,  even  though  the  general  principles  would 
almost  certainly  be  similar. 

6.  Once  it  became  possible  to  reliably  manipulate  the  BMSS  data 
so  that  geometrically  acceptable  images  using  radiance  values  could  be 
oroduced  on  a  routine  basis,  the  way  would  be  open  for  the  development 
of  spectral  analysis  procedures.  Even  though  difficulties  encountered 
in  manipulating  EMSS  data  suggested  that  time  and  funds  for  the  BMSS 
Study  would  not  be  adequate  to  attempc  an  actual  experiment  in  spectral 
analysis,  attempts  were  made  to  anticipate  difficulties  that  were  likely 
to  occur  in  that  area,  in  the  unlikely  event  that  the  data  manipulation 
problems  required  less  time  than  expected. 

7.  As  an  example,  it  was  anticipated  that  the  very  great 
differences  in  the  sensor  altitudes  and  scanning  geometries  of  the 
EMSS  and  the  BMSS  would  make  it  impractical  to  use  the  same  eorrelc'  ion 
relations  between  radiance  and  suspended  materials  concentrations  t  r 


*  To  avoid  repetition,  the  following  acronyms  will  be  used:  EMSS  =  ERTS-1 
roultispectral  scanner;  BMSS  =  Bendix  aircraft-mounted  multispectral 
scanner . 
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the  BMSS  data  as  were  used  for  the  EMSS  data.  Three  major  factors 

could  be  expected  to  contribute  to  the  differences. 

a^.  The  pixel  sizer,  of  the  two  systems  are  quite  different. 

2 

The  EMSS  pixel  has  an  area  of  about  4470  whereas, 

the  BMSS,  flown  at  an  altitude  of  about  3000  r,  has  a 

2 

pixel  covering  about  36  n  along  tne  ground  track  of 
the  flight  line.  Thus,  the  BMSS  would  "see"  much  smaller 
features  than  the  EMSS.  For  example,  spatially  complex 
features  like  eddy  systems  in  water  may  not  be  detected 
ir.  EMSS  images  but  may  be  detected  easily  in  BMSS 
Images . 

jo.  The  BMSS  scans  from  -40  to  +40  deg  with  respect  to  the 

perpendicular  from  the  scanner  to  earth,  -o  that  the  angle 
of  the  optical  axis  of  the  pixels  varies  through  a  range 
of  30  deg;  the  range  for  the  EMSS  is  less  than  10  deg. 
Thus,  the  reflectance  geometries  (i.e.  the  angular 
relation,  among  light  source,  terrain  surface,  and  sensor) 
of  the  two  systems  would  be  quite  different.  The  amount 
of  energy  reflected  by  the  same  suspended  material  c'-noen- 
trations  in  different  parts  of  the  BMSS  scan  swath  could 
well  be  different,  while  that  of  the  EMSS  could  be  treated 
as  a  constant.  There  was,  therefore,  rer..-w  table  doubt 
that  simple  correlation  relations,  as  used  for  EMSS  wata, 
would  be  useful  for  BMSS  data. 

Because  of  the  differences  in  scanning  geometries  ana 
sensor  altitudes,  the  optical  paths  taken  by  energy  cuanta 
in  the  two  system-  are  quite  different.  The  effects  of 
this  were,  difficult  to  anticipate. 

3.  Tr.  view  of  all  of  these  differences  and  uncertainties ,  there 
•  as  considerable  :c  ibt  that  an  automated  interpretation  based  on  EMSS 
.1  v  •eld  close  1-  resemble  one  made  with  EMSS  !nta. 
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Objectives  and  Scope 


Objectives 

9.  Primary  objective.  The  primary  objective  of  this  research 
effort  was  to  develop  data-handling  procedures  such  that  BMSS  digital 
data  could  be  transformed  into  radiance  values  and  used  to  produce 
images  free  of  skew  and  reflectance  geometry  distortions.  As  a  practical 
natter,  the  objective  was  to  produce  the  best  image  possible,  from  the 
BMSS  digital  data  using  information  regarding  flight  and  sensor  conditions 
recorded  at  the  beginning  of  the  NASA-prcvided  CCT. 

10.  Secondary  objective.  The  secondary  objective,  to  be  pursued 
only  if  time  and  funds  permitted,  was  to  develop  computer  software  that 
would  make  possible  the  mapping  of  the  distributions  of  spectrally 
definitive  landscape  features,  such  as  concentrations  of  suspended 
materials  in  water,  with  BMSS  data. 

Scope 

11.  Tire  research  effort  v;as  confined  to  the  use  of  existing  BMSS 
data  and  to  the  use  of  a  PDP-15  computer  for  basic  data  manipulation. 
Essential  peripheral  equipment  -was  restricted  to  an  Optronics  film 
writer  and  conventional  input-output  devices,  such  as  line  printers, 
teletype  terminals,  and  tape  and  disc  driver. 

12.  Test  areas  were  restricted  to  the  estuaries  of  the  Choptank 
River  in  Maryland  and  the  Rappahannock  River  in  Virginia.  Data  from 
only  one  flight  line  for  each  estuary  were  used  in  the  development  of 
computer  procedures. 


PART  11:  DESCRIPTION  OF  2 4 -CHANNEL  SCANNER  SYSTEM 


Mechanics,  Optics,  and  Electronics 

13.  The  NASA-operated  BMSS”*  consists  of  an  aircraft-mounted 

nul t ispectral  scanner  and  a  ground-based  data  analysis  system.  The  BMSS 
is  an  imaging  spectrometer  that  is  mounted  in  a  NC-130B  aircraft  (Figure  1). 
The  Data  Acquisition  System  (DAS)  is  a  ground-based  system  that  processes 
and  rerecords  the  spectral  data  on  computer  tapes  and  also  provides  an 
imagery  display  of  the  collected  data  on  a  viewing  cathode-ray  tube 
(Figure  2). 

14.  The  BMSS  scanner  collects  energy  in  the  spectral  range  0.34 
to  13.0  um,  which  includes  the  near-ultraviolet,  visible,  and  near- 
infraivd  regions  of  the  spectrum.  The  spectral  energy  is  separated  into 
24  channels  (or  wavelength  bands)  by  two  grating  spectrometers  sharing  a 
common  scan  mechanism.  The  scanning  mechanism  is  a  flat  mirror  that 
scans  the  earth  through  a  slit  in  the  ‘nitton  of  the  scanner  and  an  open 
door  in  the  bettor,  of  the  aircraft,  lie  mirror  rotates  counterclockwise 
with  respect  to  the  aircraft  head ing  providing  lef t-to-right  pixel 
scanning,  while  the  forward  motion  of  aircraft  provides  a  sequential 
line-by-line  scan  oi  the  earth  below.  The  scanning  mirror  views  the 
ear t :■*  through  80  deg  and  completes  a  360-deg  rotation,  during  which  it 
senses  four  calibration  sources  (Figures  3  and  -).  During  che  earth, 
scan  the  reflectance  from  each  of  700  pixels  is  processed,  and  the 
response  of  each  of  24  sensors  to  a  pixel  reflectance  is  recorded  on 

r 

magnetic  tape.  According  to  the  NASA  training  manual,  the  calibration 
sources  are  expected  to  permit  calibration  oi  the  spectral  data  to  a 
high  degree  of  accuracy.  A  precise  vitlo  of  aircraft  velocity  (V)  to 
height  (H)  ,  V/l! ,  must  b<  tai  tain  e  non  <  relation, 

i.e.  adjacent  scan  lines  just  touching  each  other  along  the  ground  track 
1  nter  if  t 1  ;<  •  wath  .  i  pi  Ld<  th i  normal  an ,  the  scan 

.mot  or  speed  i  s  ..  scablished  in  terms  of  Y/h. 
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15.  The  scan  mirror  is  set  at  a  45-deg  angle  with  respect  to  the 
aiicraft  longitudinal  axis.  The  scan  mirror  reflects  radiant  energy 
from  the  ground  and  from  the  calibration  sources  to  a  Dall-Kirkham  tele¬ 
scope.  A  roll-rate  gyro  is  used  to  maintain  the  scanner  f ield-of-view 
(FOV)  constant  regardless  of  aircraft  roll  (within  +3  deg).  The  Dall- 
Kirkham  telescope  transmits  the  reflected  energy  through  a  system  of 
mirrors  and  a  0.2286-cm  aperture  onto  a  series  of  dichroics  and  grating 
spectrometers,  which  separate  energy  of  wavelengths  less  than  2  pm  from 
those  of  2  pm  and  greater  (Figure  5).  The  aperture  is  dimensioned  to 
define  within  the  center  of  the  received  energy  field  a  FOV  of  2  miad. 
Figure  6  shows  the  path  of  the  short  wave  reflected  energy  to  reach  the 
detectors . 

16.  Scanning  in  a  lef t-to-right  rotation,  the  mirror  first  scans 
the  earth  surface,  then  the  calibration  source  elements  in  this  order: 

16  calibration  elements  from  sensor  exposure  to  the  calibration  instru¬ 
mentation  housing  (not  used),  16  from  the  UV-VIS-IR  (ultraviolet-visible 
infrared)  integrating  sphere,  16  from  the  sky  radiance  tube,  16  from  the 
high-temperature  blackbody,  and  finally  16  from  the  low-temperature 
blackbody.  Figure  7  is  a  simplified  block  diagram  of  the  complete 
system.  Figure  6  shows  the  calibration  grating  and  tining  sequence  in 
the  scan  rotation.  The  UV-VIS-IR  integrating  sphere  establishes  the 
maximum  radiance,  and  the  low-temperature  blackbody  the  minimum  railiar.ee 
for  the  scan. 

17.  The  same  detection  procedures  are  applied  to  radiance  from 
calibration  sources  as  fro:;  earth  elements;  therefore,  the  calibration 
radiance  value  fer  a  scan  and  the  recorded  pixel  radiance  from  each 
channel  are  on  the  same  basis.  Figure  9  illustrates  the  reflectance 
from  the  scan  pixels  as  the  earth  is  scanned  through  80  deg  and  the 
corresponding  reflectance  for  each  channel  in  the  calibration  sources  as 
tne  360-dog  rotation  is  completed. 

18.  The  UV-VIS-IR  sphere  (Figure  10)  is  a  16-in. -diam  aluminu"- 
sphere  painted  with  high-ref lectance  paint  and  illuminated  from  the  side 
bv  a  2C0-w  tungsten  halogen  lamp.  The  scan  mirror  views  the  illuminated 
sphere  through  a  large,  plain  optical  glass  window  and  "sees"  a 
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diffused  surface  of  supposedly  known  spectral  radiance.  The  radiance 
emitted  from  the  sphere  covers  the  range  0.3  to  2.5  urn  and  ia  calibrated 
prior  to  flight  to  a  standard  lamp  (Q3  #11)  of  known  radiance  per 
wavelength  band  (Table  1)  . 

19.  The  sunlight-matching  part  of  the  UV-VIS-IR  calibration 
contains  a  filter  wheel  with  eight  different:  neutral-density  attenuation 
screens,  which  is  located  between  the  lamp  and  the  UV-VIS-IR  sphere  and 
provides  a  means  of  changing  the  spectral  intensity  to  natch  the  sunlight 
conditions  during  a  data  run.  This  equipment  was  not  operational  during 
the  data  runs  investigated  in  this  study. 


Data  Tape  Format 


20.  The  24-channel  data  are  in  the  universal  format  as  defined  by 
NASA .  The  2 A  channels  of  data  are  digitized  into  3-bit  words  and 
recoiled  on  12  tracks  of  a  14-track  tape  recorder.  Two  consecutive 
•hansels  of  scanner  data  are  interleaved  onto  one  tape  track;  the  first 
•  referred  to  as  the  odd  channel,  and  the  second  is  referred  to  as 
•  avn.  Calibration  sources,  aircraft  information,  such  is  altitude  and 
h’c  it’.on,  and  other  housekeeping  information,  are  multiplexed  together 
on'  recorded  do  the  13th  track.  The  14th.  tra_k  contain,  the  time  code 
is .  lata  . :  •  lirect-ri  rd<  Li  anehestei  -phase-L)  code. 

Nine-track  CCT  containing  these  data  were  requested  lro::. 

record  contains  the  information  in  the  house- 


The  CCT  non 


and  to-  auxiliary  data  annotation  set  (ADAS/ASQ)  recorded  on  the 
ori  ;:unl  tapes  (figures  11  and  12).  A  program  was  written  to  present 
there  data  in  a  rare  usable  form,  (see  paragraph  31).  An  example  is 
•.haven  Ln  Table  2.  Note  the  channel  order  of  data  requested  by  the  user 
is  al  odd  ehinnuls  in  sequence  followed  by  all  even  channels  in  sequence 
(hr  i  tetel  group  in  Table  2).  This  is  the  order  of  data  on  CCT  received 
-  r->r  NASA.  Ft  j 1  owing  the  header  record  are  the  requested  data  scans. 

ntaii  7QC  :  imples ,  foil ove 3  by  BO  calibration  samples 


^  *> 


for  each  channel  (Figure  33).  Data  are  interleaved  by  bands  in  each 
sample,  two  channels  to  a  data  word. 

22.  For  seven  channels,  each  scan  on  the  CCT  contains  5460  8-bit 
bytes  of  data  (43,680  bits).  The  80  calibration  samples  are  from  the 
four  calibration  sources  (see  paragraphs  14-16)  plus  a  spare  (16  samples 
each).  Data  from  two  of  the  calibration  sources,  the  integrating  sphere 
and  the  low-temperature  blackbody,  are  used  to  establish  the  range  of 
reflectance  reception  by  each  channel  at  the  time  of  the  scan. 
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PART  111:  RESEARCH  PLAN 


Introduction 


23.  When  solar  energy  strikes  the  earth's  surface,  it  is  reflected, 
absorbed,  and  combined  with  energy  emitted  by  the  earth  (Figure  14).  As 
it  passes  through  the  atmosphere,  solar  energy  is  subjected  to  scattering 
and  absorption,  thus  limiting  the  amount  of  radiation  incident  upon  the 
earth.  The  set  of  curves  in  Figure  15  shows  the  varying  amounts  of 
absorption  at  different  wavelengths  of  the  spectrum.  The  top  curve, 
m  =  0 ,  is  the  spec'ral  curve  of  irradiance  before  it  enters  the  earth's 
atmosphere.  The  m  =  1  curve  describes  the  solar  irradiance  at  sea  level 
when  the  sun  is  directly  overhead.  The  other  curves  illustrate  the 
situation  at  various  sun  angles.  The  figure  shows  that  the  larger  the 
sun  angle,  the  greater  is  the  attenuation  of  energy.  This  results  from 
the  longer  paths  through  the  atmosphere;  the  closer  the  sun  is  to  the 
hori-on,  the  greater  is  the  amount  of  atmosphere  that  the  rays  must 


penotr  ite . 

Figure  16  shows  radiance  transmission  through  the  atmosphere, 
fhc  dips  show  the  areas  in  which  absorption  has  taken  place  due  to  water 
vapci  ,  carbon  divide,  ozone,  molecular  oxygen,  and  Rayleigh  and  aerosol 
•:  nattering. J  The  modal  wavelengths  of  each  of  the  24  channels  (or 
wave 1 eng tl  bands]  ir «  indicated  in  th<  figure  and  are  ief in<  in  Table  3. 
Vote  thac  5  of  the  24  channels  fail  within  the  visible  spectrum,  and 
that  :  1  cover  the  entire  range  of  visible  and  near-visible  wavelengths 


f0.35  to  1.1  'r..).  Ibis  study 


:s:ders  channels  4-10  which  closely 


parallel  the  run. e  of  wavelength  band  a  received  by  the  ERTS  sensors 
( Tab ic  a)  . 


25.  Clear  water  that  is  free  from  suspended  sediments  reflects 
-ost  of  the  blue-green  and  reen  portions  of  the  spectrum  and  absorbs 
the  other  visible  radiant  energy.  In  the  near-infrared  portion  of  the 
♦ectrum  (0. 8-1.1  cm),  air  energy  is  •  rbed  by  water.  A  record 

o  r  •  i  me  r or  the  earth  suri ice  in  this  portion  of  the  spectrum  can 

define  approximately  -he  land  water  :r.t e-fare  of  a  water  body.  The 


\ 


resolution  of  the  interface  is  determined  by  the  size  of  the  ground 
element  (pixel)  defined  by  the  field  of  view  of  the  scanner;  and  since 
the  total  radiance  from  a  pixel  may  be  partly  from  water  and  partly  from 
nonwater  or  from  wetted  soil,  a  clear  definition  of  the  water's  edge  is 
not  possible.  Although  absolute  definition  is  impossible,  good  approxi¬ 
mation  can  be  made  within  the  dimensions  of  two  pixels  in  length  and 
width . 

26.  Suspended  material  in  water  causes  changes  in  the  amount  of 

reflectance  and  in  the  variety  of  reflected  spectral  energy,  i.e.  the 

amount  of  energy  at  various  wavelengths.  Theoretically ,  these  changes 

are  associated  with  th *  amount  and  kind  of  suspended  material.  This 

theory'  was  supported  by  WES  investigation  of  data  obtained  from  ERTS-1 

2 

overpass  in  the  Chesapeake  Bay. 

27.  Considering  these  facts  and  theories,  a  research  plan  was 
developed  which  fell  naturally  into  two  parts  which  are  illustrated 
schematically  in  Plates  i  and  2.  Plate  ]  describes  the  plan  to  convert 
24-chanr.el  data  recorded  on  NASA  CCT  to  WES  computer  program  formats  and 
to  write,  analyze,  and  correct  geometric  distortion  when  images  are 
written  from  that  data  using  the  WES  film  writer.  Plate  2  describes  the 
plan  to  develop  the  analytical  procedures  required  to  map  spectrally 
defined  regions.  In  this  context,  a  spectrally  defined  region  (i.e.  a 
patch  of  landscape  that  displays  a  unique  combination  of  radiance  values 
in  seme  selected  set  of  wavelength  bands)  is  one  that  has  a  specified 
concentration  of  materials  suspended  in  water.  In  the  plates,  tne 

blocks  tracing  the  various  stages  are  numbered  sequentially;  the  subsequent 
discussions  are  keyed  to  the  relevant  block  numbers  in  the  plates. 

Anticipated  Problem  Areas 

28.  The  BMSS  scanning  geometry  is  quite  different  from  that  of 

the  EMSS  and,  in  consequence,  quite  different  kinds  of  geometric  distortions 
were  expected  to  occur  (Block  12,  Plate  1).  The  anticipated  distortions 
are  of  two  major  kinds. 
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£.  Scanning  geometry.  The  aircraft  was  flying  at  an  approxi¬ 
mate  altitude  of  3050  m  during  the  missions  used  in  this 
study.  Since  the  scanner  swings  through  an  arc  of  80  deg, 
the  ground  swath  covered  by  the  scanner  is  about  5118  m 
wide.  Since  the  solid  view  angle  of  the  scanner  is  fixed, 
the  pixels  at  the  edges  of  the  scan  swath  will  be  much 
larger  than  the  pixels  along  the  center  of  the  scan  swath. 

The  Optronics  film  writer,  which  is  used  at  WES  to  write 
images  from  digital  data,  writes  a  square  pixel  of  fixed 
size,  and  thus  a  direct  transcription  of  the  BMSS  data 
onto  fi'ji  with  the  film  writer  will  result  in  a  narrower 
swath  than  the  true  scanned  swath.  It  was  clear  that  this 
distortion  would  have  to  be  removed, 
b.  Aircraft  attitude.  While  aircraft  roll  was  automatically 

compensated  for  by  on-board  platform  stabilization,  aircraft 
yaw,  if  any,  could  be  expected  to  introduce  skew  into  the 
image.  Since  it  was  almost  certain  that  there  would  be 
some  wind  from  the  side,  the  aircraft  heading  would  diverge 
to  some  degree  from  the  line  of  flight.  A  procedure  would 
be  needed  for  transforming  the  BMSS  tape  data  arrays  to 
correct  for  this  effect. 

29.  The  calibration  relations  by  which  transducer  voltage  values 
(digital  form  on  (XT)  can  be  transformed  into  radiance  values 

(Block  10,  Plate  1)  in  the  BMSS  system  are  different  from  the  relations 
used  for  the  EMSS  system.  Thus,  new  calibration  transforms  would  be 
required . 

30.  One  major  problem  was  that  there  was  r.c  possibility  of  deter¬ 
mining  the  actual  suspended  materials  concentrations  in  the  waters  of 
the  test  areas  during  the  time  of  the  aircraft  overflights.  This  would 
have  required  teams  in  hosts  collecting  water  samples,  and  this  was  far 
too  costly  to  be  practical.  As  a  result,  there  was  no  possibility  of 
establishing  an  empirical  relation  between  suspended  materials  concentration 

and  radiance  (Rlock  18,  Place  2)  as  had  been  done  previously  for  the 

2 

studies  using  the  EMSS.  The  best  that  could  be  hoped  for  was  a 
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relation  that  described  relative  differences  in  suspended  materials 

concentrations,  based  on  the  fact  chat  low  concentrations  of  suspended 

2 

materials  appear  to  produce  direct  linear  correlations  with  radiance. 
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PART  IV:  DEVELOPMENT  OF  DATA  HANDLING  PROCEDURES 

Data  Conversion  Software 

Blocks  1  and  2,  Plate  I 

31.  With  the  help  of  Reference  6,  a  computer  program  was  written 
that  decoded  the  BMSS  CCT  and  reassembled  the  contained  data  into  forms 
more  readily  usable  by  WES  computer  hardware  and  available  software. 

Several  modules  were  involved.  For  example,  the  NASA  tapes  contain 
housekeeping  data  (Figure  11) ,  aircraft  location  and  attitude  data 
(Figure  12),  radiant  energy  data  from  terrain,  and  calibration  data. 

These  data  were  all  reformatted  so  that  they  could  be  printed  out  in  a 
more  readable  form,  as  illustrated  in  Table  2. 

32.  Another  module  was  designed  to  transform  the  scanner  transducer 
data  (i.e.  the  output  of  the  sensors  that  measure  the  amount  of  energy 
received  by  each  BMSS  wavelength  band)  into  the  format  required  by  the 
WES  Optronics  film  writer.  This  format  is  illustrated  in  Figure  17. 

The  general  principle  is  that  a  new  tape  is  devoted  to  four  wavelength 
bands  or  dummies  so  that,  each  wavelength  band  can  be  processed  independently 
with  existing  WES  software. 


Test  Site  Selection 


Blocks  3  and  4,  Place  1 

33.  In  the  prior  work  with  ERTS-1  data,  maps  were  produced  of 

suspended  materials  concentrations  in  several  estuaries  leading  into 
2 

Chesapeake  Bay.  Of  the  several  available,  the  Rappahannock  and  Choptank 
estuaries  were  chosen.  The  P.a  p  pah  anno  ck  River  estuary  was  chosen  because 
it  has  sources  of  relatively  uniform  suspended  materials  and  a  small 
tidal  range.  For  contrast,  the  Choptank  River  estuary  was  chosen  because 
its  sources  of  suspended  materials  are  quite  complex,  and  its  _idal 
range  is  somewhat  greater  than  that  of  the  Rappahannock.  The  locations 
of  the  selected  study  areas  are  shown  in  Figure  18.  EMSS  tapes  covering 
these  areas  were  ordered  from  NASA. 
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Data  Transformation 


Blocks  5,  6,  and  7,  Plate  1 

34.  Once  the  programs  to  read  and  reformat  the  NASA  BMSS  tapes 
had  been  written,  the  NASA  tapes  covering  the  BMSS  records  for  the  two 
test  areas  were  reformatted  (Block  1)  into  WES  Optronics  film  writer 
format  (Figure  17).  At  this  stage,  the  contained  data  have  not  been 
transformed  in  any  way.  The  data  relevant,  to  scanner  pixels  are  still 
raw  digital  values  defining  transducer  output  voltages,  and  the  array  is 
a  simple  and  rigid  raster  array.  If  used  to  drive  the  Optronics  fiSm 
writer,  the  resulting  image  would  include  all  of  the  distortions  intro¬ 
duced  by  scanning  geometry,  aircraft  yaw,  and  so  on. 

35.  In  addition,  all  housekeeping  data  (Block  6),  calibration 
data  (Block  7),  and  so  on,  were  extracted  and  assembled  into  tables  as 
illustrated  in  fable  2.  Among  the  critical  information  are  those,  items 
describing  flight  parameters.  Of  particular  importance  are  the  tine, 
location,  radar  altitude  (which  is  a  measure  of  absolute  height  above 
the  terrain),  true  heading,  drift,  and  date.  The  significance  of  these 
parameters  will  be  discussed  later. 

Water-Land  Discriminat ion 


Blocks  3  and  9 ,  Plate  1 

36.  Because  the  secondary  objective  (see  paragraph  10)  involves 

the  distributions  of  suspended  materials  in  water,  it  was  obviously 

advantageous  to  delete  .ill  nenwater  areas  fren  the  reformatted  BMSS 

tapes,  since  this  would  result  in  a  significant  saving  of  processing 

time.  To  achieve  this,  the  water-land  discrimination  program  developed 

2 

tor  use  of  EMSS  data  was  used  without  modification.  When  used  with 
EMSS  data,  the  discrimination  between  water  and  land  is  made  with  the 


rear-infrared  band  (ERTS-1  channel  7;  0.8-1. 1  un)  ,  since  ail  energy  in 
tuts  wavelength  band  is  nearly  totally  absorbed  by  water,  while  it 
reflects  strongly  from  vegetation,  soil,  and  other  "land"  materials. 
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37.  The  properties  of  F.RTS-1  channel  7  are  most  closely  approximated 
by  those  of  BMSS  channel  10  (0.981-1.045  yin.  Table  3),  and  this 
wavelength  band  was  therefore  used  to  define  the  pixels  relating  to  open 
water  surfaces.  This  information  was  then  used  as  a  "digital  mask"  (see 
procedure  described  in  Reference  2)  to  edit  out  all  "land"  areas  in  all 

of  the  wavelength  bands  used  in  this  study  (Block  8)  (3MSS  channels  4- 

10)  . 

38.  Despite  the  fact  that  the  imagery  resulting  from  writing  the 
raw  tapes  is  badly  distorted,  it  is  helpful  to  write  an  image  of  each 
wavelength  band  (Block  9),  simply  to  confirm  that  all  decoding  and 
reformatting  procedures  have  operated  correctly  (Figure  19).  Even  an 
uncritical  comparison  of  the  image  with  a  map  of  the  corresponding  area 
(Figure  20)  ..-ill  reveal  that  the  geometry  is  badly  distorted. 

Voltage-to-Radiance  Transformation 


Blocks  10  and  11,  Plate  1 

39.  Using  the  reformatted  data  and  the  calibration  data  contained 
on  the  EMSS  tapes,  a  program  was  written  that  transformed  the  values 
recorded  by  the  BMSS  into  radiance  values.  Since  calibration  data 
varied  considerably,  no  attempt  wan  made  to  develop  a  conversion  equa¬ 
tion  in  which  constants  were  applicable  to  the  entire  data  set.  Instead, 
a  conversion  based  on  the  data  from  calibration  sources  from  each  scan 
was  derived  and  applied  to  that  scan  only. 
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=  radiance  at  earths  surface  of  i  pixel  in  channel  j 
=  recorded  digital  value  of  reflectance  from  i^  pixel  in 
channel  j 

=  average  of  16  samples  from  the  1 ov-temperature  blackbody  for 
channel  j  for  a  scan 

*=  average  of  16  samples  from  the  integrating  sphere  for 
channel  j  for  a  scan 

=  percent  reflectance  with  which  the  integrating  sphere 
compares  to  the  standard  lamp  for  channel  j 
=  spectral  radiance  evaluation  of  the  standard  lamp  per 
micrometer  for  channel  j 
=  channel  width  of  j  channel 
=  atmospheric  transmittance  factor 
-  sun  angle  from  nadir  (perpendicular  from  sensor 
to  earth) 


For  each  data  collection  run.  (o^)  ((RSL^)(B^))  (I/O  (cos  SA)  for  each 
channel  is  constant. 


•  Image  Geometry  Correction 

Blocks  12,  13,  14 ,  and  15,  Plate  _1 

40.  Two  major  kinds  of  geometric  distortions  occur  in  the  BMSS 
data;  namely,  that  caused  by  aircraft  crabbing  or  yaw  (i.e.  the  angular 
relation  between  aircraft  heading  and  the  ground  track  of  the  aircraft) 
and  that  caused  by  the  scanning  geometry  (see  paragraph  28a) . 

41.  The  procedure  used  to  correct  for  aircraft  crabbing  is 
illustrated  in  Figure  21.  Line  A ' B '  in  Figure  21a  indicates  the  actual 
scan  relation  to  the  center  line  (i.e.  the  true  course  of  the  ground 
track  of  the  aircraft)  of  the  scene.  The  orthogonal  grid  array  of  data 
on  the  BMSS  tape  places  the  scan  line  at  90  deg  to  the  scene  center  line 
(line  AB  in  Figure  21a  and  b),  which  causes  a  misplacement  of  features 
as  indicated  in  Figure  21c.  The  displacement  at  each  end  of  the 

scan  can  be  computed  from  the  difference  in  the  true  heading  and 
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true  course. 


By  staggering  blocks  of  pixels  across  the  scar,  in  a 
•diagonal  array  down  scans  (Figure  21d  and  e)  ,  features  aie  adjusted  to 
their  normal  relation  to  each  other  and  to  the  center  line  of  the 
scene.  The  geometries  illustrated  by  Figure  21n-e  were  reduced  to  a 
computer  program  that  rearranges  the  pixels  in  an  orthogonal  grid  such 
that,  when  written  with  the  film  writer,  the  skew  created  by  aircraft 
crabbing  is  eliminated.  An  example  of  this  type  of  correction  is  given 
in  Figure  22.  Note  that  alter  the  correction,  the  ends  of  the  views  are 
not  at  right  angles  to  the  center  line. 

42.  A  quite  different  procedure  is  required  to  compensate  for  the 
variation  in  pixel  width  ns  the  instantaneous  viewing  angle  is  removed 
from  ti.°  center  J  ine  (Figure  2’i).  At  nadir  the  instantaneous  FOV 
(2  nrad)  is  6.1  m  (20  ft)  wide  when  the  aircraft  is  flying  at  an  altitude 
of  30  4.6  a  (10,000  ft).  The  24-channel  Forth  observations  Aircraft 
Program  flight  over  the  Rappahannock  on  22  April  1973  vas  flown  at 
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O 


~  (id, 550  ft),  which  defines  the 


oixe 


width  at  nadir  as  6.6/ 


ft),  bach  .successive  \  •  .-.el  is  .:jder  as  th*.  scanning  angle  from 


a  r.  an 


no  : nereis es 


pixel  viewed  at  the  scan  extreme  is 


of  -0  deg  the.  center  line  and  is  11.31 


-.  vice, 


■  re  '  . :  the  at  fix  ■  taj  lata  tc  Lnd  ti  the  patia] 

ri  itions  between  th<  :arth  pixels  from  which  data  ;.,ore  drawn,  some 

Lve  proced  >plied  tii  rrays  was  necessaiy.  To 

s:  r<.-tch  the  scan  >rrcy  a.  data,  pixel  radiance  values  ware  repeated  at 

hese  points  !  •  ring  the  difference 

in  distances  from  the  nadir  point  when  computed  using  th.o  tangent  of  the 
scanning  angle  and  the  altitude  of  the  sensor  and  when  computed  using 
the  number  of  pixels  from  nadir  and  the  width  of  the  pixel  at  nadir. 

..lie:'  the  difference  in  these  distances  exceeded  the  width  cf  the  pixel 

.  '  radii  aJ  le  wa  repeated  expanding  the  number  of  pixels 

in  the  scan.  This  expanded  scan  pushes  the  scan  extremes  to  more  nearly 
coincide  with  the  earth  location  sensed  by  the  extreme  instantaneous 
r,.v'.  In  Figure  2J,  pixel  75  has  been  repeated.  The  farther,  away  from 
the  renter  :  in.:,  th?  more  often  pixels  are  repeated.  In  the  Rappahannock. 
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digital  nap,  70  pixels  were 
of  this  situation  were  also 
of  a  scene  before  and  after 
Figure  24.  Note  that  after 
than  before  correction.  It 
to  write  an  example  of  the 
the  programs  have  operated 


added  to  each  half  of  each  scan.  The  geometries 
reduced  to  a  computer  program.  An  example 
correction  for  scanning  geometry  iB  given  in 
correction  the  scene  is  significantly  wider 
is  usually  helpful,  as  illustrated  in  Figure  24, 
scene  in  each  wavelength  band  to  verify  that 
properly. 
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PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusions 


objective 


•'*3.  The  primary  objective  was  achieved.  A  set  of  computer  programs 
is  now  available  that  transform.  the  NASA-provided  BMSS  CCT  into  an 
image  formed  of  radiance  values  and  free  of  both  skew  and  scanning 
geometry  errors.  All  programs  are  run  «  n  a  ILU-13  computer,  and  images 


are  produced  directly  with  an  a  routes  film  v 
L'lSS  tapus  can  be  used  effectively  by  agencies 
computer  facilities  to  create  pictures  of  o r co¬ 
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The  Bendix  24-channel  sensor  has  been  abandoned  In  favor  of  the  Bendlx 
11-c.hannel  sensor  because  of  the  nore  consistent  performance  of  components 
less  sensor  "noise",  and  better  control  of  the  internal  flight  environment 
The  channels  investigated  in  this  study  are  included  in  the  Bendix  11- 
channel  sensor  equipment  with  the  same  channel  wavelength  ranges  and 
same  type  receptors.  It  is  recommended  that  this  study  be  repeated  with 
existing  data  collected  by  the  Bendix  11-channel  sensor  over  water 
bodies . 
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Table 


Calibration  Data  for  Standard  Lamp  (QB  ’( 1 1 ) 


Wave  length 

Spectral 

Radiano 

um 

w/sr/ pm/ cm- 
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X 

1 
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r—i 
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6.48 

X 

io"4 
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1.22 

X 
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X 

I 

o 
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2.87 

X 

IQ'3 
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X 

10~3 
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(7) 
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X 
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X 

I 
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X 
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(ID 
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X 

1 

O 
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1.3 

(24) 
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X 
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5.33 
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1.9 
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X 

io" 3 
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(27) 

1.28 
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(26) 

2.5 
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5.41 

X 
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Table  2 


Information  extracted  From  Rappahannock  River  Tapis 


Record  Sire  156 
Nc.  of  Elements  700 
Channel  No.  5 
Channel  No.  7 
Channel  No.  9 
Channel  No.  4  * 

Channel  No.  6 
Channel  No.  8 
Channel  No.  10  _ 

Scan  StatuG  255 
Starting  Scan  Line  66351 
Time  16:46:15.3 
Latitude  N37  deg  37.3  rdn 
Longitude  W  76  deg  28.5  min 
Radar  Altitude  10950 
Barometric  Altitude  36800 
True  Heading  Deg.  285.5 
Drift  Deg.  L  Pnan  Course  1.4 
Roil  Rt .  Wing  D  Deg.  0.7 
Pitch  Nose  D  Deg.  0.2 
Ground  Speed  271  kts 
Late  4-22-73 
Mission  230 
Site  244 
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Table  3 


BMSS  Channel  Characteristic  Bands  Compared  to  EMSS  Bands 
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Figure  1.  Airborne  raultispectral  scanning  systeri'* 
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Figure  C round -based  data  anal  vs  is  system 
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Figure  j.  Schematic  diagram  of  collecting  opti 
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Figure  7.  Simplified  basic  multispectral  ser.sor  function  diagram 


V \C0VP 
k  NCO0 C* 
'  •  T  $ 


IUIM  rCAL< 
GATES 


BACKG  ../JND 
GATES  (NOT  USEOi 

CAL  R(AD CVGS 


JAG  ppr 


SPATE 


f  .  V'VlS 


ca'fs 


J  S-vV  Rif 
'fS  <  - 


JL 


JL 


;•  .  *  \  '.n(!'.\0  i 

i 


( — 


SOT  <JSI  Z" 


A  ‘i  I  AO 
;  4  .s***  >\(>S 


J  L 


-o+  ^  ^  »■*  *•  •  3*} 


Figure  9.  Composite  spectral  input 
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Figure  11.  BMSS  housekeeping  frames 
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Figure  12.  ADAS/ASQ-90  packing  in  hoi  sel  eepir.g  voi  ds 
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Figure  14.  Interaction  of  radiant  energy  with  matter 


OX(m)  WATTS  M  PER  MILLIMICROMETER 


Figure  15.  Solar  spectral  irradiance  curves  at  sea  level  with 
varying  optical  air  masses^ 
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STEP  1  -  Extract  data  in  more  usable  form  from  NASA  CCT. 
Extract  channels  5,  7,  S,  and  10  on  one  tape 
and  4,  6,  8,  ar.d  10  on  another.  Include 
appropriate  calibration  data  at  end  of  each  scan. 
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STEP  3  -  a_.  Write  two  intermediate  classification  tapes 
using  EKTS  program. 

b.  Using  a  special  24-channel  program,  merge 
the  two  intermediate  tapes  to  produce  one 
classification  tape. 
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V..tS  procedure,  showing  re 


Figure  17. 
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STEP  4  -  Apply  24-channel  program  to  correct  for  crabbing 
(essentially  a  skewing  technique  in  the  direction 
of  line  of  "flight).  A  new  tape  is  produced. 


STEP  S  -  Apply  24-channel  program  to  expand  scan  image 
producing  a  new  tape. 
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STEP  2  -  a. 

b. 


c. 


Rewrite  tape  data  in  format  for  ERTS  programs 
on  hand,  masking  land  if  desired. 

Apply  24-channcl  algorithm  to  convert  digital 
tape  data  to  radiance  at  earth  surface  as  new 
tapes  are  vrritten  (if  desired). 

Produce  desired  digital  pictures  from  any 
channel  data  by  coding  the  digital  data 
(0^255)  to  256  shades  of  gray  on  the  Optronics 
film  writer. 


c t_  (see  FlGoAG  2b) 
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STEP  6  -  Write  and  develop  film  Interpreting  digital 
codes  into  shades  of  gray. 


Optronics  film  writer  from  2k -channel  data 


a.  FLIGHT  LINE  3  -  LOWER  RAPPAHANNOCK  RIVER.  VIRGINIA 


Figure  18.  Flight  line,  locations  (Mission  218  and  Mission  230) 
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Figure  19.  Film  written  from  digital  CCT  data  (0-255),  Rappahannock 
River,  Mission  230,  22  April  1073 
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Figure  21.  Data  correction  procedure  for  aircraft  crabbing 
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Flight  altitude 


(a) 


Earth  surface  represented  by  tape  pixels 


Figure  23.  Correction  procedure  for  registering  scan  extremes 
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Primary  Objective:  Plan  to  Develop  Data  Handling  Procecares 
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PLATE  1 


Secondary  Objective.  Plan  to  Pevelon  Corputer  Programs  to  Map  Suspended 

Materials  Concentrations  with  SMSS  Data 
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PLATE  2 


In  eeccrdfjsse  with  ER  70-2-3,  persgrapb  6c(l)(b), 
data d  15  Fabruery  1573»  e  facsimile  catalog  card 
in  Library  of  Coojresa  format  Is  reproduced  below. 
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